Explosion-implosion duality in the Bose-Einstein condensation 



Pijush K. Ghosh* 
Department of Physics, Ochanomizu University, 
2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, Japan. 



m 
o 
o 

c 
a 
>—» 



> 
o 

On 
O 



o 
o 



We show an explosion-implosion duality in the one and 
two dimensional Bose-Einstein condensation with or without 
a particular time-dependent harmonic trap. The result is in- 
dependent of the strength and the attractive or the repulsive 
nature of the self-interaction of the condensate. This implies 
that the implosion in a particular atomic species has a dual 
description in terms of the explosion in the same or another 
atomic species and the vice versa. The result is applicable 
without any modifications to non-relativistic theories that are 
invariant under the SL(2, R) reparameterization of the time- 
coordinate. 

PACS numbers: 03.75.Fi, 05.45.Yv, 03.65. Ge 



Duality plays an important role in science. Models 
with a dual description are also abundant in different 
branches of physics. The duality in the two dimen- 
sional Ising model [1], the electromagnetic duality in su- 
persymmetric gauge theories [2] , the explosion- implosion 
duality in astrophysics [3,4] are a few celebrated exam- 
ples among them. The purpose of this letter is to show 
an explosion-implosion duality in the Bose-Einstein con- 
densation (BEC) within the framework of the Gross- 
Pitaevskii equation(GPE). 

The action describing the two dimensional BEC with- 
out the harmonic trap is invariant under a SL(2, R) repa- 
rameterization of the time coordinate, leading to a dy- 
namical 0(2, 1) symmetry of the model [5]. The same is 
true for the one dimensional BEC with repulsive interac- 
tion [5]. One special case of the SL(2,R) reparameter- 
ization of the time coordinate corresponds to the inver- 
sion, r — ► t = -i. The field and the space-coordinates 
also transform simultaneously [5] to keep the action in- 
variant. This particular symmetry, which we refer to as 
the duality-symmetry, is the central to the discussions of 
this letter. We show that the duality-symmetry leads to 
explosion-implosion duality in the GPE without the trap. 

The presence of the harmonic trap is essential for many 
experimental observations in BEC. Thus, it is natural 
to study if the duality-symmetry discussed above can 
be preserved in the BEC even in presence of the har- 
monic trap. Unfortunately, the introduction of a time- 
dependent harmonic trap breaks both the scale invari- 
ance and the invariance under a translation in time. 
This leads to breaking of the dynamical 0(2, 1) symme- 
try. We, nevertheless, show that the duality-symmetry is 
present for a very particular choice of the time-dependent 
harmonic trap. The time-dependent frequency of this 
particular trap, which varies inverse-squarely with time, 



is independent of the space dimensionality and the same 
for both the one( repulsive case only ) and the two di- 
mensional BEC. 

We further show that the invariance under the duality- 
symmetry leads to explosion-implosion duality in the 
one( repulsive case ) and two dimensional GPE with the 
particular time-dependent trap. In particular, following 
the method of Ref. [5], we first obtain the exact time- 
dependence of the width of the wave-packet in the phys- 
ical problem described in terms of r. We also obtain 
the same in the dual problem described in terms of vari- 
ables obtained from the physical model by the duality- 
transformation. The width in the physical problem is 
zero initially. The width spreads gradually with the in- 
crease of the time and finally diverges at the asymptotic 
infinity. On the other hand, in its dual description, the 
width diverges at t = — oo. As t increases, the width first 
takes a finite, but, large value and then decreases with 
the passage of time. The width finally vanishes at t — 0. 
Thus, the explosion in the physical problem has a dual 
description in terms of implosion in the dual problem, 
leading to an explosion-implosion duality in the BEC. 
We remark here that, unlike in the case of wave collpase 
at a finite time [6], the explosion/implosion occurs only 
at the final value of the allowed ranges of the time coordi- 
nate r/t. This difference is very crucial in understanding 
our results. 

The consequences of such an explosion-implosion dual- 
ity is the following. The self-interaction in the GPE can 
be either attractive or repulsive depending on whether 
the s-wave scattering length of the atomic species is neg- 
ative or positive, respectively [7]. Further, using the Fes- 
chbach resonance method [8] , the magnitude and the sign 
of the self-interaction between atoms of a certain species 
can be tuned to any desired value : large or small, re- 
pulsive or attractive. If the self-interaction is repulsive 
and sufficiently strong so as to overcome the effect of 
the confining potential, the condensate is expected to 
explode. On the other hand, we expect the condensate 
to implode in the strongly attractive regime of the self- 
interaction, if the zero-point energy due to the kinetic 
energy term can no longer balance the attraction. How- 
ever, contrary to this known behavior of the condensate, 
we show that the exact explosion-implosion duality is in- 
dependent of the strength and the attractive or the re- 
pulsive nature of the self-interaction. Thus, it is possible 
to have both explosion and implosion for any fixed value 
of the scattering length of an atomic species, provided 
the condensate is prepared satisfying the required initial 



1 



conditions described below. The reason behind such an 
interesting, novel and apparently counter-intuitive pre- 
diction is that there is no scale in the problem to prevent 
an initially growing/collapsing condensate from further 
growth/collapse. The explosion in an atomic species can 
also be described in terms of an implosion at a differ- 
ent value of the scattering length of the same or another 
atomic species and the vice-versa. 

Consider the following non-relativistic Lagrangian in 
arbitrary d + 1 dimensions, 

c = ird T ^-^-\ wi 2 - ff mnr). a) 

Am 

The coupling constant g has the inverse-mass dimen- 
sion in the natural units with c = h = 1. The po- 
tential V is real and does not depend on any dimen- 
sional coupling constant. This implies a scale- invariance 
in the theory. We demand the invariance of the action 
A = J drd d vC under the following time-dependent trans- 
formations [9-11], 



r -> r h = f(t) ^ r , t -» t = t(r), f(t) 
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with the scale-factor r given by, 
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Particular choices of r(t) = t + (3, a 2 t, and 
spond to time translation, dilation and special conformal 
transformation (SCT). The Noether charges correspond- 
ing to these symmetry transformations close under an 
0(2,1) algebra [5]. 

The action A is invariant under the transformation (2) 
and (3) for g = [10]. For g ^ 0, the invariance of A 
under (2) and (3) solely depends on the form of V. We 
choose, 



v(\ v |) = | V 
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(4) 



This choice of V gives the critical non-linear Schrodingcr 
equation(NLSE) in d+1 dimensions [6]. The Lagrangian 
(1) with the above choice of the interaction has an 0(2, 1) 
symmetry for arbitrary d. Note that d = 1 gives a sex- 
tic interaction that is relevant for one dimensional BEC 
with the repulsive interaction [12]. The BEC in 2 + 1 di- 
mensions is described by a quartic nonlinear interaction 
[7]. Note that for d = 2, we indeed have V = | V t ■ This 
produces a cubic NLSE which is also directly relevant for 
two dimensional optics [6,11]. As an aside, we remark 
here that for d = 3, the equivalent formulation of (1) and 
(4) in terms of hydrodynamic variables [13] is directly 
relevant for the supernova explosion or implosion in laser 
induced plasma [3]. 



The Lagrangian (1) with the potential given by (4) is 
invariant under the transformations (2) and (3). Con- 
sider a very particular case of Eq. (3), 
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This implies the following symmetry-transformations for 
the field ip and the coordinate r, 
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ip(T,r) ^ ip h (t,r h ) = - exp[i—r ]tp{T,r), (6) 
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which is known as the lens transformation in the context 
of the critical NLSE [6,11]. The parameter (3 is real and 
arbitrary. Consequently, the theory at a time r > is 
mapped to a theory at a time t < and the vice versa. 
In particular, the physical problem at r = is mapped 
to the dual problem at t = — oo. Similarly, the physi- 
cal problem at r = oo is mapped to the dual problem 
at t = 0. The critical value separating this two regime 
is obviously given by r = or t = 0. Thus, the phys- 
ical problem is bounded from below in time, while the 
dual problem is bounded from above in time. The above 
mapping also involves a time-dependent scaling for the 
space-coordinate. We choose the scale- factor 4 = — ^ 
to be positive-definite and follow the convention for the 
time-coordinates and the parameter (3 as given below, 



< t < oo, -oo < t < 0, (3 < 0. 
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Note that (3 has been chosen to take negative values only. 
If one prefers to choose a positive value for (3, the ranges 
of the time-coordinates in (7) should also be interchanged 
simultaneously. We remark that the positivity of the 
scale-factor is essential for the consistency of our analysis. 
For example, the densities p = ip*^ an d Ph = V'/lVVt are 
related to each other through the following equation, 



Ph(t,T h ) 
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This is consistent for odd d only if the scale-factor is 
positive-definite. Further, the width of the wave-packet 
( see below ) in arbitrary d becomes negative, unless we 
demand the positivity of the scale-factor. 

Let us now introduce a moment / and its dual X as, 



J(r) = f J d d v r 2 P , lit) = ^J d d r h 
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The moment / and its dual X can be interpreted as the 
expectation value of the square of the radius of the con- 
densate. Using the transformation (6), it is easy to see 
that the moment / and its dual X are related to each 
other by the equation, 
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Following Ref. [5] and Eq. (10), 7(r) and X(t) are uni- 
versally given by, 
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where a and b are arbitrary constants. The constant of 
motion Q = IH — (5^7) 2 is related to the Casimir opera- 
tor of the underlying 0(2, 1) symmetry [5]. Note that all 
the information on the system under consideration ( like 
the strength and the attractive or the repulsive nature of 
the interaction, the space-time dimension on which the 
problem is considered, the nonlincarity etc. ) is con- 
tained in the expression of 7 and X in the constant of 
motion Q only, through the Hamiltonian H . Thus, the 
dynamics can be described in terms of the same set of 
initial conditions for any value of g: positive or negative, 
large or small. It is possible to have both explosion and 
implosion for a particular value of <?, if the condensate is 
prepared following the respective initial conditions, which 
can be extracted from the exact expressions of 7 and X. 
One might naively think this as unphysical. However, 
note that both the kinteic energy term and the interac- 
tion term scales in the same way. So, there is no scale 
in the problem to prevent an initially growing/collapsing 
condensate from further growth/collapse. 

The criteria for the collapse of the condensate at a fi- 
nite and real time r* is Q < 0. In particular, the moment 
7 vanishes at a finite time r* , 



(a 2 b 2 + Q) 
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which is real if Q < 0. Note that we have the freedom of 
making r* either positive or negative by choosing appro- 
priate values for the integration constants a and b. The 
moment 7 is semi-positive definite by definition. Thus, 
the exact expression for Q = 7i7— (A_^7) 2 implies that the 
condensate collapses for any initial condition if if < 0. 
On the other hand, if H > 0, the condition for the col- 



lapse is given by, 
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explosion-implosion duality is forbidden for Q < 0, since 
there can not have any explosion either in the physical 
or in the dual theory. 

We now consider the case Q > for which the 
explosion-implosion duality can be realized. A positive 
Q necessarily implies a positive H, H > j{\^) 2 - This 
lower bound on H can be equivalcntly written as a con- 
straint on the initial profile of the condensate, 
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where 77 is the dual Hamiltonian. Using the exact time- 
dependence of 7, X and an alternative expression for 
H(H), 77 = \§?z{'H = 50), it is easy to verify that 
such a criteria is indeed satisfied. Finally, whether a con- 
densate satisfying the criteria (13) will explode or im- 
plode depends on whether the expectation value of the 
square of the radius of the condensate is finite or infinite, 
respectively, at the initial time. 

The condensate in the physical problem is of finite ex- 
tent at r = with 7(0) = a 2 . The condensate swells 
gradually as r increases and explodes at r = oo leading 
to a divergence in 7. In the dual description, the conden- 
sate is of infinite extent at t = — oo with a divergent X. As 
t increases, the condensate is extended over a very large, 
but, finite area. It then shrinks with the advancement of 
time, finally reaching to a constant value for the moment 
X at t = 0, 2/(0) = f3 2 (b 2 + . Thus, the evolution of 

the condensate in the physical problem is different from 
that of its dual description. In particular, the explosion 
in the physical problem has a dual description in terms 
of implosion. Note that we have the freedom of choos- 
ing a 2 > y/Q\ p \,b 2 = (af3)- 2 (a 4 - Q/3 2 ), such that, 
1(0) = 7(0) = a 2 . This allows the condensates to have 
the same radius, both at the initial time in the physical 
problem and at the final time in the dual description. 

A comment is in order at this point. The lens trans- 
formation (6) has the property of converting a collapsing 
solution to a stationary solution known as the Townes 
soliton [6,11]. The Hamiltonian vanishes, when evaluated 
for the Townes soliton [6] , implying that the constant of 
motion Q is necessarily negative. So, there is no contra- 
diction between the known results on the Townes solition 
and the results of the present letter on the explosion- 
implosion duality which is strictly valid for Q > 0. 

We now show that it is possible to choose a particular 
time-dependent harmonic trap maintaining the duality- 
symmetry. Consider the following Lagrangian, 

C h = C-^mu(T)r 2 \i>\ 2 1 W (r) - c 2 ( 7 r - a)' 2 , (14) 

where we have made a very specific choice for the fre- 
quency of the time-dependent trap and ojo is an arbitrary 
constant. The action Ah — J d,Td d r£h is invariant under 
the transformation (2) with the r(i) given by, 



r(t) = 



at + 13 
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(15) 
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Note that the time-translation, the dilatation and the 
special conformal transformation can not be obtained 
as special cases from the above relation. However, the 
transformation responsible for the duality-symmetry is 
contained in (15). This is obtained by putting a = 
and, hence, 7 = — 4, 



W (r) 



(16) 
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For this choice of the time-dependent frequency, the ac- 
tion Ah is invariant under the transformations given by 
the Eqs. (5) and (6). Thus, although the action Ah is 
not invariant under the individual time-translation or the 
SCT, it is indeed invariant under the duality- symmetry. 

Using the transformations (2) and (3), the action A can 
be mapped [5] to the action Ah with the time-dependent 
frequency given by (16), if we choose, 



r(t) = 



1 
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where c\ and ci are arbitrary constants. The parameter 
i] is determined in terms of (3 and cjq, 
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For | ujq | > 2pr, V becomes a complex number. More- 
over, at | coo | = 2pp T W becomes divergent. Thus, we 
restrict ujo to take values within the range, < | wo | < 
2py. Consequently, r\ is restricted to as, \ < rj < 1. Note 
that Ah is invariant under the transformations (5) and 
(6) for arbitrary loq. The restriction on luq is valid only 
if we want to relate A to Ah through a time-dependent 
coordinate transformations. The BEC without the trap 
can be obtained in the limit — > or r\ — ► 1. The trans- 
formation (17) reduces to a symmetry transformation for 
i] = 1 and keeps A invariant. 

The moment / for the BEC with the harmonic trap 
can be interpreted as the square of the width of the wave- 
packet [14]. Following [5], the moment / and its dual X 
for the trapped BEC are determined as, 
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We have chosen ac 2 = ±VQ(1 — 2r/)~ 1 and ^- = ±\/Q 
in deriving Eq. (19). Recall that the coordinate t is 
allowed to take values in the negative-axis only and Q > 
0. 

The width I h vanishes at r = in the physical prob- 
lem. This is because the confinement is too strong to al- 
low any non- vanishing width in this limit. The strength 
of the confinement decreases with the increase of r and 
the width starts spreading over the accessible area. This 
gradually leads to a divergence in the width at the asymp- 
totic infinity, where the confinement is so weak that it 
can not hold the condensate any longer. This signifies an 
explosion. In the dual description of this problem, the 
width l h diverges at t = — oo. As t increases, the width 
first becomes large, but, finite and then decreases grad- 
ually. This leads to a collapse of the condensate with 
vanishing width at t = 0. This signifies an implosion. 



Thus, the explosion in the physical problem has the dual 
description in terms of implosion. 

In conclusion, we have shown the explosion-implosion 
duality in the one (repulsive case) and two dimensional 
BEC with or without the particular time-dependent har- 
monic trap. This implies that the implosion in a par- 
ticular atomic species has the dual description in terms 
of the explosion in the same or another atomic species 
and the vice versa. With the recent advancement of 
the technology related to the BEC [8,15], such a dual- 
ity may be realized in the laboratory in near future. Fi- 
nally, we would like to mention that the method and the 
results of this letter are applicable without any modifica- 
tions to non-relativistic theories that are invariant under 
the SL(2, R) reparameterization of the time coordinate. 
Physically interesting models with such an invariance are 
abundant [16-18]. It is surprising to note how the under- 
lying dynamical 0(2, 1) symmetry gives a universal result 
for models ranging from astrophysics to models as diverse 
as BEC. 
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